Leading Infrared Logarithms from Unitarity, Analyticity and Crossing 
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We derive non-linear recursion equations for the leading infrared logarithms in massless non- 
renormalizable effective field theories. The derivation is based solely on the requirements of the 
unitarity, analyticity and crossing symmetry of the amplitudes. That emphasizes the general na- 
ture of the corresponding equations. The derived equations allow one to compute leading infrared 
logarithms to essentially unlimited loop order without performing a loop calculation. For the im- 
plementation of the recursion equation one needs to calculate tree diagrams only. The application 
of the equation is demonstrated on several examples of effective field theories in four and higher 
space-time dimensions. 
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I. INTRODUCTION 

Effective Field Theories (EFTs) are non- 
renormalizable field theories, which allovif the inves- 
tigation of the infrared (low-energy) behaviour of 
various physical systems (see recent review Q). The 
standard tool for the studies of the asymptotic be- 
haviour of renormalizahle field theories is the method of 
renormalization group equations (RGEs). In the case 
of EFTs the method of RGEs must be modified as the 
number of counterterms increases rapidly with the loop 
order. 

A possibility of the systematic construction of RGEs 
for non-renormalizable^ quantum field theories was 
demonstrated in Ref. [2|. In particular, it was shown 
that the series of the leading logarithms (LLs) can be 
obtained by calculation of one loop diagrams. However, 
the solution of the RGEs derived in Ref. [2] requires the 
calculation of non-trivial one-loop diagrams, the number 
of which is rapidly increases with the loop order. There- 
fore, the implementation of this method in practice is not 
an easy task. The method of Ref. 2} has been applied 
in Ref. @ for the calculation of the five-loop LLs for the 
pion mass in the massive 0{N + 1)/0{N) sigma model. 
In Ref. Q the authors, using dispersive methods, cal- 
culated the three-loop LLs to tttt scattering in massless 
Chiral Perturbation Theory (ChPT). 

Recently, a completely different method for the calcula- 
tion of LLs in a wide class of non- renormalizahle massless 
field theories was developed in Refs. [1, The non- 
linear recursion equations derived in Refs. @ allow 
one to obtain the LLs contributions without performing 
non-trivial loop calculations at each loop order. 

In the present paper we show that the non-linear re- 
cursion equation for LLs is a consequence of analyticity, 
unitarity, and crossing symmetry of the 5-matrix. The 
requirement of analyticity, unitarity, and crossing sym- 
metry for massless particles allows to obtain LLs coef- 
ficients without calculation of loop integrals. [We note 
that these requirements in conformal field theories lead 
to very powerful bootstrap equations, see The form 



of the corresponding equations is such that one can easily 
generalize them in various ways - for arbitrary space-time 
dimension, different symmetry groups, different spins of 
the fields etc. 

To make the presentation of our method transparent 
we shall restrict ourselves to a massless scalar EFT, which 
is described by the following generic action: 
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(1) 



where the expansion of the interaction part of the ac- 
tion V{<j),d(l)) starts with four fields 0, the corresponding 
"0*" part contains 2fc derivatives. The index a (it can 
be a multi-index as well) corresponds to a possible inter- 
nal symmetry of theory ((TJ . EFTs of type ((Ij have the 
property that the anomalous dimension of fields is zero 
at the leading order. The theories with "0'^" interaction 
can be considered in similar framework, this a more com- 
plicated case will be considered elsewhere. We stress that 
the absence of mass is crucial for our discussion below. 



II. GENERAL METHOD 



We consider the 4-particle scattering amplitude 0" -|- 
0'' — ?> (/)'^ -I- (j)'^ in four dimensional space (generalization 
to an arbitrary number of dimensions is given in Sec- 
tion IVip . The LLs for other physical quantities, such 
as form factors , (generalized) parton distributions [1] , 
etc., are related to the LLs of the 4-point amplitude. 

The 2 — )■ 2 amplitude can be decomposed in the irre- 
ducible representations of an internal symmetry group of 
the EFT H]): 



A'''"'''{s,t,u) 
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A\s,t,u). 



(2) 



Here Pf"'^ 



is a projector on an invariant subspace cor- 
responding to the irreducible representation / of an in- 
ternal symmetry group. According to the Wigner-Eckart 
theorem the corresponding projectors can be obtained as 



2 



the convolution of two Clebsch- Gordon coefScients, see 
e.g.i. 

The renormahzabihty of a theory depends on the di- 
mension of the couphng constants that enter action ([1} . 
Let us consider the case of the action ([T]) with one cou- 
phng constant of dimension (4— _D— 2fc) {2k is the number 
of derivatives in the interaction part of the Lagrangian, 
D is the dimension of space-time). We denote the corre- 
sponding couphng constant as the dimension of the 
constant _F is fc-|-r'/2 — 2 (the constant F corresponds to 



where S — (^^^^pyz is a dimensionless expansion parame- 
ter, is a couphng in the Lagrangian. The scale pa- 
rameter /i is arbitrary in the LL approximation, since its 
change influences the next-to-leading logs only. ©(NLL) 
stands for terms with next-to- leading logarithms (NLL). 
The series of LLs for the amplitude ^ takes into account 
that the scattering amplitude has both right and left cuts 
in the complex plane of the Mandelstam variable s. 

Our aim is to derive the recursion relations for the 
coefficients in front of LLs in the low-energy expansion 
of the partial wave amplitudes: 

n— 1 

Obviously one has: 

n-l 
i=0 

where the index n corresponds to the loop order plus one, 
and the angular momentum I is restricted to I < kn. I 
denotes the irreducible representation of the correspond- 
ing internal symmetry group. 

In order to derive the recursive equations for the LL 
coefficients we shall employ the unitarity, analyticity 
and crossing symmetry of the scattering amplitude. 

The crossing symmetry relates the amplitudes with 
the interchanged momenta to the amplitudes of different 
symmetry group representations (the summation over the 
repeated "representation indices" is always assumed): 

A'is,t,u) = Cli'A'\s,u,t), (6) 
A'is,t,u) = Cii'A'\u,t,s). 

The crossing matrices C are defined as 

CII' Tjabcd Tjbacd ^ /^II' rjabcd rybdac ^ (n\ 
tU — ^I' ^! ^SM — ^1 \') 



the pion decay constant in usual ChPT) . The case of sev- 
eral couplings we consider in Section |Vl At fc = 2 — Z5/2 
theory ([T]) is renormalizable, and at /c > 2 — £'/2 it is a 
non-renormalizable low-energy EFT. 

We shall work with the partial wave amplitudes (see all 
definitions in Appendix A) as they depend only on one 
energy variable s. For D — A (see generalization for an 
arbitrary space-time dimension in Section IVl)) the partial 
wave amplitude in LL approximation can be represented 
in the most general form as follows: 



(3) 

I 

where di is the dimension of the irreducible representa- 
tion / (see the explicit expressions for the crossing ma- 
trices below). 

The partial wave scattering amplitudes possess right 
and left cuts in the complex plane of the Mandelstam 
variable s. The right cut discontinuity is fixed by the 
unitary relation: 

Imif(s) = |tf(s)P + ©(Inelastic part). (8) 

The unitarity relation is valid for the physical region 
s > 0. Under the inelastic part we understand the part 
with more than two particles in an intermediate state. 
Obviously the LLs coefficients do not depend on this part 
of the expression. Indeed, by cutting more than two lines 
in any diagram the power of logarithms decreases and 
hence such cuts influence only NLLs. We can restrict 
ourselves to the elastic unitarity relation, since only LL 
coefficients are considered. 

Substituting the LL expression ^ for t\ (s) into the 
unitarity relation (jS]) and collecting the coefficients in 
front of S^L"~'^ {L = In (^^/s)), we find the relation 

n — 2 — 1 ^ n— 1 

{n-i- l)a^;- = ^ ^Wn~rd- (9) 

Thus, the unitarity allows us to relate a and lo coeffi- 
cients. However, corresponding relations © do not allow 
us to obtain the closed form equation for the coefficients 

LO. 

We need some additional information about the left 
cut. This information is provided by the dispersion rela- 
tions (analyticity) and the crossing symmetry ([6]). One 
can show that if only two-particle cuts are taken into ac- 
count, the following relation connects discontinuities on 
the left and right cuts (see derivation in the Appendix A) 



J 



oo n— 1 2 2 

«; w - f i: ^ E -a (^) 1""-- (^) + °(NLL), 

n=l j=0 



3 



l'=0 " 



This relation for the tttt scattering is weU-known, it is 
a consequence of the Roy equation 10], it was used in 
Ref. [4] to calculate three-loop LLs in ChPT. Substituting 
Eq. ^ into Eq. (fTO|) and collecting the terms in front of 
the LLs we obtain the following relation in addition to 

®: 

n — i — 1 kn n—1 

E - E 2FfT E-^-"-M'(-i)'^''"i;«'(ii) 

i=o ;'=o 2=1 

where the matrices fiJJ are defined as follows 

I 



[Here we introduce the matrix Cst which is defined as a 
product of the crossing matrices d?]) Cst = CsuCtuCsu]- 
The recursion relation (|12l) for the LL coefficients should 
be supplemented by initial conditions, i.e. by the values 
of dX n — 1 and correspondingly ^ = 0, 1. The corre- 
sponding values can be obtained by the tree level calcu- 
lation of the partial wave amplitudes t\ (s) (see Eq. @) 
using action ([T]) of an EFT under consideration. 

The recursion equation p2p is the main result of the 
present paper. This equation allows the calculation of 
the leading infrared logarithms to an essentially unlim- 
ited order. Furthermore, this method presents a puissant 
tool for the study of the general structure of the infrared 
logarithms. For the derivation of Eq. (fT2|) we used only 
the unitarity, analyticity, and crossing symmetry of the 
amplitude. This fact emphasizes a general nature of the 
non-linear recursion equation (|12p . The corresponding 
equation can be derived for many physical problems de- 

I 



where cr^ ^ i^^ _ Y,a=i 't^'^'P"- The part does not 

contribute to the 4-particle amplitude due to the absence 



I 

The x (n + 1) matrices fiJJ perform the crossing 

transformation for the partial waves. Main properties 
and explicit expressions for the matrices fi" are given in 
the Appendix B. 



Summing relations © and (ITT]) and taking into ac- 
count t -O- M crossing symmetry we arrive at the follow- 
ing closed non-linear recursive relation for LL coefficients 



(12) 

I 

scribed by a non-renormalizable effective low-energy La- 
grangian, e.g. theory of critical phenomena, low-energy 
quantum gravity, theory of magnetics, etc. 

In the next Sections we apply our method to several 
EFTs. Also we give the generalization of the method on 
an arbitrary dimension and on the theories with mixed 
renormalizable and non-renormalizable interactions. 

III. 4D MASSLESS 0{N + 1)/0{N) cr-MODEL. 

We have several reasons to consider the 0{N -|- 
l}/0(iV) cr-modcl. First, this model was considered in 
[5| where the recursion relations for LLs were derived by 
a completely different method. Second, this model (at 

= 3) is equivalent to the Weinberg Lagrangian [ll| of 
ChPT. The Lagrangian of the 0{N + l)IO{N) cr-model 
has the form 



(13) 

I 

of masses. The cr-model (fT5]) belongs to the here consid- 



n — 1 kn ^ 

^ E E E ^ ('^"'^" + ^"^in + ciii-ir'^Z 



J , ,J 



=1 r=o 



21' 



c = -{d.ad'^a + d.ro^r) - ^d^rd^r - ^(</>'^r)a2(//) + 0(</.«), 
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ered class of EFTs ^ with D = 4 and k=l. 

The theory possesses 0(iV) symmetry, therefore 
the expansion in invariant tensors ([2|) goes over the ir- 
reducible representations of the 0{N) group. Namely, 
over the three representations of dimensions dj — 



1, 



jV(Af-l) {N+2){N 



— I . We enumerate these represen- 



tations by index / = 0,1,2; such naming corresponds 
to the value of the isospin in the s-channel for the case 
iV = 3. 

The corresponding projection operators on these rep- 
resentations have the following form: 



^abcd 



ab scd 



N 



Pi 



abed 



ac s:bd 



gadgbc 



The straightforward calculation with help of Eq. ([7]) gives 
the crossing matrices 







o\ 




Ctu — 





-1 


c { 




^0 


1/ 






Substituting the crossing matrices into Eq. (1121) we obtain 
explicit form of recursion relations for the LL coefficients. 
The initial conditions for the recursion can be obtained 



J 



abed 



2 V 



^bd ^ gadgbc 



ab s:cd 



N 



by the trivial tree-level calculation of the scattering am- 
plitude with the help of the Lagrangian (jTS]). The result 
is (jjIq^ = N — \, ojIq"^ = —1, and wff ^ = 1, aU other 
coefficients w(j are zero. 

At first glance, the recursion relation (IT^ is different 
from that derived in Ref. Q. The reason is that the LL 
coefficients w^j with fixed 0{N) quantum numbers are 
not independent and can be expressed in terms of the 
universal LL coefficients ujni which were used in Ref. [H . 
The corresponding relations have the following form: 



l'=0 



l'=0 



l'=Q 



(14) 



The origin of such relations is very simple - the 0{N + quantum numbers are interrelated. 

1)/0{N) a model is a theory with one coupling constant The equation on uj is obtained from (IT2|) by the sub- 

and hence the channels with different internal symmetry stitution (jl4[) . it reads 

I 



m=l i=o 1=0 
even even 



with the initial condition wio = l,wii = .The cocffi- 

j_ j-y(mA)(n—m,l) ■ i 

cients B- '^ ' are given by: 



(m,i)(n-m,l) 



B 



(m,i){p,l) _ 



2j + l 



i^ + i-m E 



min[p,m] (^ikr^lkr^k] 



fe=0 



2fc + 1 



(15) 



(16) 



Eq. coincides exactly with the equation obtained Ref. [5| . In the present paper the equation was obtained 

by a completely different, more complicated method in 
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FIG. 1: The lowest partial wave LL coefficients in the 
0(4)/0(3) sigma model up to the 137th loop order. 



from a simple symmetry consideration. The detailed dis- 
cussion of the properties of the recursion equation (jlSp 
(and hence of Eq. (fT^ 'l and its particular solutions can 
be found in Ref. [5[. For convenience of the reader we 
present in Appendix C [Tables I-III] the values of the LL 
coefficients uj^i up to the 4-loop order. In Fig. [T] we show 
the values for the lowest partial waves LL coefficients in 
the 0(4)/0(3) sigma model up to the 137th loop order. 



LL coefficients are restricted only to the lowest partial 
wave ^ = 0. Using this simplification and introducing the 
following notation o;,^ = oj^^q we can rewrite the recur- 
sion relation for the renormalizable theory p7)) in a 
simpler form: 



,/=o 



1 



1 ^ 

i=l 



N- 



2N 



(19) 



N -2 



2N 



.1=1 



.1=2 



n — 


1 ^ 




i 


1 




n — 


1 ^ 




i 


N - 


- 1 



N 



1 

N 



(20) 
(21) 



This system of equations can be easily solved with the 
following result: 



/=0 



N 



.1=2 



-{N + 2)[- 



N + % 



.1=1 



0. 



IV. 4D 0{N) SYMMETRIC 0* THEORY 

As a next example we consider a renormalizable field 
theory. For simplicity we consider the same 0{N) sym- 
metry as in the previous case. A possible renormalizable 
Lagrangian of type dU) is: 



(17) 



Due to the same 0{N) symmetry as for the cr-model 
the crossing matrices are the same. The differences to the 
previous case are the number of derivatives in the inter- 
action vertex, i.e. A; = in Eqs. (|4ll2p and the coupling 
constant F is dimensionless now. 

The low-energy expansion of the scattering amplitude 
in the renormalizable theory (jl7p goes over increasing 
power of LLs (note that S = l/(47ri^)^ at fc = in 
Eq. Q): 



Substituting this solution into the LL expansion of the 
scattering amplitudes (jl8|) . we can perform the summa- 
tion of the LLs with the result: 



TT iV + 2 



1 



,1=2 



2 (47rF)2 

TT 2 

2 (47rF)2 



bi 



T hi 



(23) 



hi 
(47r_F)- 



■In 



All other amplitudes are zero in the LL approximation. 
In Eq. (1231) the constant bi = -^^^j^ is the well known 
result for the one-loop beta function of the 0{N) sym- 
metric 0^ theory ([TT]), see e.g. Ref. [l^ . 

The logs enter Eq. (j23p in the combination in which 
one can easily recognize the running coupling constant 
(denoted as in this paper) for the 0{N) symmetric 

(/)^ theory. It is not surprising. Indeed, for a generic 
renormalizable theory (e.g. k = 0, D — 4) the recursion 
equation is reduced to the following form: 



2 ^ 21 

n—1 



1 (47rF)2 



.ln"-i(^). (18) 



The LL coefficients w^^; satisfy the recursion relation 
with k — Q. The initial conditions are obtained by cal- 
culation of the scattering amplitude at the tree level in 



theory pT|) - w 



/=o 

10 



-(iV-f 2) 



,1=2 



-2 with all other 



n — 1 



i=l J 



(24) 



where the matrix B^"^ is independent of the index n and 
it is expressed in terms of the crossing matrices: 



B'' = -{l + Cst + Csu) 



1 



ij 



(25) 



Eq. p^ is simplified considerably at fc = as it con- 
tains only the matrix il" = (5"^(5'°. This implies that the 



The recursion equation can be reduced to a simple 
system of differential equations. For that we introduce a 



6 



generating function for the LL coefHcients as follows 
uj^{t) = J2'^=i^n t^~^- Then the recursion relation ([M)) 
is reduced to the following differential equation: 



dard RGEs in the case of a renormalizable field theory. 
Eq. can be considered also as the generalization of 
RGEs for the case of non-renormalizable EFTs. 



dt 



(26) 



with the initial conditions that can be obtained by the 
calculation of the tree diagrams, e.g. in 0{N) symmet- 
ric (j)^ theory oj^ {0) = {-N - 2, 0, -2}. Eq. ® has the 
form typical for the RGEs for the running coupling con- 
stant. In this way, we demonstrate that our recursion 
relation (IT2]) for LL coefficients is reduced to the stan- 



V. THEORY WITH RENORMALIZABLE AND 
NON-RENORMALIZABLE INTERACTIONS 

In this section we consider the theory with mixed 
renormalizable and non-renormalizable interactions. To 
be specific we consider a Lagrangian which is the sum of 
the 0{N ^ l)/0(iV) cr-model and the 0{N) symmetric 
(/)^ theory: 



J 



1 



1 



r 



+ (27) 



where = ~ '^a=i 4'°'4'°'- In this case the expansion 
is performed over two parameters: a dimensionless cou- 
pling A and dimensional We note that at iV = 3 
the theory (|27p corresponds to the famous two-flavour 
chirally symmetric Weinberg Lagrangian plus the chiral 
symmetry breaking term proportional to the constant A. 
That chiral symmetry breaking term is proportional to 
the constant ^3 in the notations of Gasser and Leutwyler 
[isj . However, here we force this term to be of the same 
order as the Weinberg Lagrangian. 

The series in infrared LLs for the amplitude in the the- 
ory (P?]) can be presented in the form of double expansion 
in the parameters and A: 



Here S = 



n.m— 



21 + 1 



2 

yn^^n+m-l (^l^Y (28) 



i^^^PY and A — A/(47r)^ are dimensionless ex- 
pansion parameters. The LL coefhcicnts w^^; now have 
an additional index m, which shows the order in the 
constant A. Also, for convenience, we define Wqq; = 0. 
Obviously, w^q; coincide with the LL coefficients for the 
0{N + \)/0{N) (T-model considered in Section Hill and 
correspondingly the coefficients Wg^^ are equal to the LL 
coefficients for the 0{N) symmetric 4)^ theory (see Sec- 
tion |IV|. 



Imposing the conditions of the unitarity, analyticity 
and crossing symmetry on the LL amplitude (I28p and 
repeating steps described in Section [Hi we arrive at the 
following recursion equation for LL coefficients (^nmi- 



in 



n.m n 



ir s:IJ 



^ 2 



+ ci,;:(-i)'+' a 



r 



^ijV ^n—i,m—j,l' 

2FT1 



(29) 



As usual, the initial conditions for this recursion equa- 
tion can be obtained by the tree level calculation of the 
scattering amplitude with the Lagrangian (l?7l) . Actually, 
the initial conditions are the combination of the corre- 
sponding initial conditions discussed in the previous two 
sections. 



An important feature of the recursion equation (j29|) is 
that its dependence on the index m is trivial. Eq. ((^^ 
with fixed indices n and / has the generic form (1^^ . 
Therefore, the corresponding equation, in principle, can 



be solved for fixed indices n and /. Introducing the gen- 
erating function: 



(30) 



we can partially sum up the LL expansion of the ampli- 
tude lli): 



'/(«)-fi:".'4^K^))5^i""-'(?)(3i) 
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This expansion is similar to the LL series for the theory 
with a single coupling constant (U), the only difference 
is that the LL coefficients w^^j are functions of A and the 
log. Expressions for the functions are rather com- 

plicated, therefore we give only a couple of (non-trivial) 
examples for = 3, corresponding to ChPT. We con- 
sider the case of n = 2, that corresponds to the one- loop 
chiral log which is dressed by any number of loops with a 
interaction. The result of simple calculations gives: 



.=o.,^ _ 950 160 10 

100 , 10 

A i H 

221-'^ ' 39' 



'iT\t) 



^/.(t) + |^/3W, (32) 

100 



with the functions fi{t) given by: 



m) 



1 



(33) 



where hi = 11/2 is the coefficient of the one- loop /3- 
function of the 0(3) symmetric (j)^ theory. The "anoma- 
lous dimensions" 7^ have the following values: 71 — 
9/11,72 = 40/33,73 = 10/33. Using the result ^ one 
can see that the inclusion of the (j)^ interaction to the 
massless chiral Lagrangian (0(4)/0(3) cr-model) leads 
to decreasing of the coefficients in front of the one-loop 
chiral log (except the J = 1 amplitude for which the one 
loop chiral log is zero). Possible physics implication of 
this will be considered elsewhere. 



VI. EFFECTIVE FIELD THEORIES IN 
ARBITRARY DIMENSION 

The derivation of the recursion equation (jl2p was done 
for D ~ A. It is straightforward to generalize it for an ar- 



bitrary even dimension D > 4 [l5| . We restrict ourselves 
to the even dimensions in order to avoid the appearance 
of power-like singularities in the amplitudes. The ana- 
lytic structure of the amplitude in an odd dimensional 
theory should be investigated separately. 

We consider the Lagrangian of type ([1]) with the in- 
teraction part containing 2k derivatives. For D > 4 all 
theories of type ([T|) are non-renormalizable. Simple di- 
mensional analysis shows that an {n — l)-loop diagram 
in a theory ([1]) scales with low external momenta as 
^ (p2-|fen+(n-i)(_D-4)/2^ rp]^-g power bchaviour is accom- 
panied by logs with the maximal power of {n — 1). The 
low-energy expansion of the partial wave amplitude t\ (s) 
[details of partial wave decomposition in D dimensions 
can be found in Appendix A] in LL approximation has 
the following form: 



•,D-4 



32ns- 



E 

n=l 



21 + D -3 



.(34) 



We introduce here the dimensionless expansion parame- 
ter: 



(47r)^2 



D 



F2 



(35) 



Obviously, the expansion (p4| is reduced to Eq. ^ for 
D = 4. 

Now, with the help of Appendix A the reader can easily 
repeat the derivations presented in the Section |lT] and 
after simple calculations arrive at the following recursion 
equation for the LL coefficients a;^^; for a £)-dimensional 
theory: 



n-i kn+{n-l) 



1 ^ 

1=1 



E 

i'=0 



- 

2 



fcn+(n-l)- 



fcn+(n-l)- 



21' + D-3 



.(36) 



The crossing matrices fi" (D) are straightforward gener- 
alization of fi" for an arbitrary dimension: 

(^)"--^(f^)^i:«::(-)c,^w. 

1=0 

Here C\ ^ (z) are Gegenbauer polynomials, which form 
a basis for the partial wave decomposition of the ampli- 
tude in D dimensions. Details on the crossing matrices 
fijj (D) can be found in Appendix B. 



Various examples of multi-dimensional theories will be 
considered elsewhere. Here, as an example, we give only 
results for the LL coefficients in the 6D 0{N) symmetric 
(/)^ theory, which corresponds to fc = 0, 1? = 6 in Eqs. ([Ml 
I36p . The LL expansion for the partial wave amplitude in 
this theory has the following form: 

n— 1 ^ 
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The four-loop results for the coefBcients w^j are presented 
in Tables IV- VI of Appendix C. 

The recursion equation (|36|) for the 0{N) symmetric 
(p* theory in D dimensions can be solved in the large N 
limit. In this limit the amplitude is dominated by the 
S-wave and by the singlet / = "isospin" component. 
Details of the large N limit for the recursion equations 
can be found in Ref. ^] . The result for the LL coefficients 



-N 



N 



2{D-3) 



(38) 



With help of Eq. (|34p we can perform the large- TV sum- 
mation of LLs in the D-dimensional 0{N) symmetric (f>* 
model: 



327rr (^) (D - 3)F2 
1 



(39) 



1 



(47r) T 2^-3 (I5-3)r( )F2 



In ^ 



VII. CONCLUSIONS 

Using the requirements of the unitarity, analyticity and 
crossing symmetry of the scattering amplitude, we derive 
recursion equations for the coefficients in front of leading 
infrared logs in massless effective field theories. The cor- 
responding equations are given by ()12[) for AD theories 
and by (|36p for an arbitrary even dimension. To imple- 
ment these equations one needs to perform a calculation 
of tree diagrams only. One needs such calculation to find 
initial conditions for the recursion equations (|12l36p . 

In Section [ni] we demonstrate that our recursion equa- 
tions for LLs (| 121) are equivalent to the recursion equation 
derived in Ref. (5|1 by a more complicated method. The 
method of Ref. [5| required non-trivial all-order analy- 
sis of the structure of possible counter-terms as well as 
rather complicated loop calculations. In the present pa- 
per we found a much simpler and more general way to 
derive the corresponding recursion equations. 

In Section IIVI we show that our recursion equations 
(|12l36p are reduced to usual RGEs for the case of a renor- 
malizable field theory. Therefore, one can consider our 



recursion equation for LLs as a generalization of RGEs 
for the case of a non-renormalizable field theory. 

Eq. (fT2|) in the case of a renormalizable field theory 
is reduced to a simple first order differential equation 
(RGE). For a general case it can be reduced to a Ham- 
merstein integral equation [l3 | . For that type of integral 
equations one can prove that the solution does exist and 
is unique. A little is known about explicit solutions of 
the Hammerstein integral equations, although many ap- 
proximate methods have been developed. One of them 
is the reduction of the integral equation to the recur- 
sion equation of the type (fT2|) . Indeed we checked on 
many examples of EFTs that the corresponding recur- 
sion equation is very effective and it allows one to obtain 
LL coefficients to essentially unlimited loop order (e.g. it 
take a couple of minutes on a PC to compute the 99-loop 
LL coefficients for massless ChPT). 

The derivation of our main equations is based on gen- 
eral properties of a quantum field theory: the unitarity, 
analyticity and the crossing symmetry. The specific form 
of the theory enters the equation only through the cross- 
ing matrices (type of internal symmetry) and the initial 
conditions for the recursion equation. That shows that 
the recursion equation for leading infrared logs has a gen- 
eral nature and can be easily written for any massless 
EFT with fields of various spins (e.g. gravity or theory 
with spinor fields) and in an arbitrary dimension. 
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Appendix A: Analytical properties of partial waves. 

In this appendix we give a summary of definitions and 
relations for partial waves in arbitrary dimension D. 

The partial waves in D {D — 1 spatial plus 1 time) 
dimensions are defined as 



A\s,t) = 647r^ 



21 + D -2,^[^] 



Ci {cos e)tj{s) 



(Al) 



T{l + D~3) Jo 



0-3 2t 

d9 sm°'^eA\s, COS e)C I ^ {cosO), cos6i = H-— , 
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where Cl'{z) are Gegenbauer polynomials. Note that 

Ci^^{z) = Pi{z), hence the expansion (jAl|) for Z? = 4 is 
reduced to the usual partial wave expansion in Legendre 
polynomials. The unitarity of the ^-matrix leads to the 
following relation for the partial scattering amplitudes: 



Imtf(s) = TT^S 



©(Inelastic part) 



|tf(s)p (A2) 



This unitarity relation allows one to obtain the discon- 
tinuity of the amplitude on the right cut (s > 0) in the 



complex plane of the Mandelstam variable s. 

The discontinuity on the left cut can be obtained with 
the help of dispersion relations. If one takes into account 
only cuts related to a two-particles intermediate state 
(that is enough for the LL approximation) the analytical 
properties of the amplitude are quite simple. The am- 
plitude has the s-channel cut from 4m^ to -|-oo, and the 
u-channel cut from to —oo. [We switch on the masses 
for a moment in order to avoid the problems with coa- 
lescing of branch points]. The usual dispersion relation 
at fixed t with no subtractions can be written in the form, 
see e.g. [loj . 



J 



ds'i 



KlI' 



C. 



II' 

su 



— s — iO s' — 4rn^ + t + s — iO 
I 



\lmA''{s\t), 



(A3) 



where the matrix Csu is the crossing matrix ([6|) . In prin- 
ciple, we have to make subtractions in this dispersion 
relation, but the subtractions do not influence the imag- 
inary part of amplitude, and therefore for our considera- 
tion we can drop them. The discontinuity of the ampli- 
tude (IA3p on the left cut (s < 0) receives the contribution 



from the second term of Eq. (|A3[) only. Computing the 
corresponding discontinuity on the left cut and perform- 
ing the partial wave decomposition (jAl[) we find the fol- 
lowing relation between discontinuities on the right and 
left cuts: 



imtf(.) = Y^di' 



jj, 2"-\2l' + D -i) 



;'=o 



s — Ani^ 



(s — Am?Y 



J 



^(^) 



r(;-i-D-3) TT 

As^Arn^ - s ~ s') 



(A4) 



c; 



/ s + 2s' - Am^ \ _ / 2s + s' - Am 



Am? 



C„ 



Am^ 



This relation for D = A for the case of the tttt scattering obtain to the expression pO|) . 

amplitude was derived in Ref. [3| as a consequence of the In the massless limit it is more convenient to rewrite 
Roy equation. Taking D — A and the limit — > we Eq. \AA\ in the form 



i'=0 

1 



T{l + D-'d) TT 



(A5) 



r 



ii' 



Appendix B: Crossing matrices Q,'. 

In D dimensions the definition of the crossing matrix 
in the partial wave space is the following: 



Using the orthogonality relations for the Gegenbauer 
polynomials, we can write the crossing matrix as the fol- 



(^)"c,-(f±f)^i:«;'(«)c,;*=w- m 



l'=0 
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lowing integral: 



2i' + D-3 2^-H'\ r^r^ ri 



T{1' + £> - 3) 



This integral can be computed in terms of hypergeomet- 
ric function in the Saalschutz form: 



(B2) 



J 



f2« {D) = 



(-l)''+"(2r + £'-3) n\ T{l + D-2,) 
T{n + l' + D- 2) (n - A\ 



1-F3 



-l,l + D-2,,-l' -n- D-3,l' -n 

D-4 D-2 
2 ' 2 



-n, — n - 



(B3) 



This representation of the crossing matrices fi" {D) is 
very convenient for the numerical calculations. 

The main properties of the crossing matrices (D) 
are the following: 



Appendix C: Tables for LL coefficients 



n 

^(-i)^nJ^(D)nf = (-i)'+''n«'(£>), 



(B4) 



0Fr (D + n - 2) 



In this Appendix we present Tables for LL coefficients 
for the AD 0{N + 1)/0{N) a-model (Tables I-IIl) and 
for the &D 0{N) symmetric (j)^ model (Tables IV- VI). 
The empty entries in the Tables correspond to zeros. 



TABLE I: Table of / = LL coefficients for the 4D cr-model, (jj'„T° • (iV - 1)"^ 



n\l 





2 


4 


1 


1 






2 


N 1 

2 


5 

18 




3 


(UN 1 59 
4 144 ' 144 


13JV 1 13 
144 ' 48 




4 


JV-^ 631JV^ 1 46279Af 13309 
8 2700 1 194400 194400 


173^^^ 4313JV 1 5333 
21fi0 38880 1 .3Rf!R0 


AT^ 49JV 1 8 
200 5400 "T R75 


5 


JV" 136N^ 1 2498743Ar^ 
16 675 ' 7776000 

3083771iV 1 619889 


1417Ar^ 1 481367JV^ 

40320 ' 3628800 
727373iV , 1071107 


N^ 1 9787W^ 

280 756000 
449681]V J 81007 


11664000 ' 4665600 


4082400 ' 6531840 


27216000 ' 5443200 



11 



TABLE II: Table of 7 = 1 LL coefficients for the 4D tr-model, w^^^ 



n\l 


1 


3 


5 


1 
± 


1 






2 


JV 1 3 

9 9 






3 


9JV^ 37JV I 49 
40 SO SO 


JV^ 37iV 1 49 
40 720 + 7-;0 




4 


JV-^ 1 493]V'^ 41791JV , 8543 
10 + 1200 64800 1 15960 


JV* 1 79iV'^ 7019iV 1 233 
40 900 64800 5592 




5 


5JV* 52859Ar^ , 18963533^^^ 
112 302400 ' 64432000 
9B8BB87Af , 2686037 


5JV* 98743Ar'^ , 4018577]V'^ 
288 1666200 ' 34992000 
3612281Af 1 292127 


Af* 20753JV-' 
2016 10886400 
1 363091JV2 17849Af , 101 


27216000 ' 10886400 


34992000 ' 4665600 


97977600 4898880 ' 40824 



TABLE III: Table of 7 = 2 LL coefficients for the 47) cr-model, w^f^ 



n\l 





2 


4 


1 


-1 






2 


AT 1 1 

^lJ_2 


iV 5 

_s_ia 




3 


1 3JV 59 
8 + 16 144 


iV^ _,_ 47JV 13 
8 ^^ 48 




4 


AT^ 857iv2 1 21131JV , 13309 
20 10800 1 194400 ' 194400 


JV^ 449JV^ 1 68711JV 5333 
14 2160 1 575160 38880 


JV^* 41JV^ 1 407iV 8 
580 3600 1 55500 675 


5 


iV* 1 1727iv3 3323209JV2 
48 ~'~ 25920 23328000 
1 1492651JV 619889 


25]V* 1 112891]V3 4774289iV^ 
672 ' 725760 16329600 
1 612299JV 1071107 


JV* 1 21797]V3 1747919JV^ 
224 1209600 54432000 
1 282487JV 81007 




"T" 11664000 4665600 


2041200 6631840 


9072000 6443200 



TABLE IV: Table of 7 = LL coefficients for the 673 ^•*-model, coif" 



n\l 





2 


4 


1 


-(iV + 2) 






2 


i(iV + 2)(iV-l) 






3 


-^(iV + 2f(iV-l) 






4 


(iV + 2f(iV-l)(^ + ^) 


-^{N + 2fiN-l) 




5 


-{N + 2f{N -1) 

n'^ I 223]V 1 1 
V 1296 151200 ' 226800-' 


-(N + 2f(N-l) 

23JV 1 101 \ 
^ 388800 ' 583200 


~{N + 2f{N - I) 

N , 23 N 
^ 680400 ' 5103000 > 



TABLE V: Table of 7 = 1 LL coefficients for the 67) 0*-model, w^f ^ 



n\l 


1 


3 


5 


1 









2 








3 








4 


(iV + 2f(^-^) 


(A^ + 2) (, r,-,,„ 58350 ) 




5 


(N + 2]^(^^'^ 1 ^ ) 


(N + 2Yi \ N 127 N 




\i\ T ^) V 27216 ' 4320 136080/ 


\iy -r ^; V 46360 ' 43200 2268000/ 
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n\l 



TABLE VI: Table of / = 2 LL coefficients for the 6D 



^ ^ >^ S4n ^ 7.7m 22a 



AT 



' 4536 



i^iv -r :;>if;n :^4?inn ^ 4Sfin 



(^ + 2)(-^+ 453600 
17JV I 1 \ 
56700 ^ 113400 



(iV + 2)(- 



11664 
I N I 101 N 
~'~ 97200 ~'~ 291600 



(7V + 2)( 

AT 



+ 



408240 
23 
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